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m ■ Abstract 



A conjecture is presented for the thermal one-point function of boundary operators 
in integrable boundary quantum field theories in terms of form factors. It is expected to 
have applications in studying boundary critical phenomena and boundary fiows, which are 
relevant in the context of condensed matter and string theory. The conjectured formula 
^ \ is verified by a low-temperature expansion developed using finite size techniques, which 

can also be used to evaluate higher point functions both in the bulk and on the boundary. 

\^ ■ 1 Introduction 

a\ 

Cp ■ The aim of the present work is to calculate the thermal one-point function of local boundary 

operators in integrable boundary quantum field theories. Such a theory can be specified with 
^ ! a Euclidean action of the form 

00 
O 



A=rdT(l dx /:($", 9^$", a^$°)+£B($"(x = 0),9^$°(x = 0))') (1.1) 

where the field variables are denoted $". The bulk equations of motion follow from the Euler- 
^ ■ Lagrange equations specified by £, while the boundary condition is obtained by varying JZb] 



the possible choices for the action are restricted by requiring integrability [T]. 

For a finite temperature T the Euclidean time r must be compactified to a volume 

T 

Consider a local operator O inserted at the boundary x = as shown in figure II. 1[ The 
quantity of interest is the thermal average 



^^ _ Tr (e-^^O) 
Tr (e-«^ 



(Of = ',^.^J (1-2) 



where H is the Hamiltonian corresponding to the action p.ip and the trace is taken on the 
space of states allowed by the boundary condition. 
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Figure 1.1: The finite temperature boundary quantum field theory with a local boundary 
insertion O 

The main motivation to study finite temperature correlators of boundary operators comes 
from boundary renormalization group flows, where the most useful quantity characterizing 
the space of the flows is the Affleck-Ludwig (7-function or boundary entropy [2]- The original 
setting where this function was introduced already made use of finite temperature. Further- 
more, as shown by Friedan and Konechny |3], the variation of this function along the flow 
can be computed via a sum rule that is expressed in terms of finite temperature boundary 
two-point functions. The present paper can be considered as a step towards constructing such 
correlators from field theory data. In addition, quantities like the thermal average (|1.2|) may 
have direct physical relevance to condensed matter systems. 

Our goal is to express the thermal average in terms of matrix elements (form factors) of the 
operator O. Therefore in section 2 the boundary form factor bootstrap is presented, slightly 
extended from its original formulation in ^ to include theories with more than one particle 
species. In section 3 we formulate a conjecture for the thermal average (11.2(1 based on the 
earlier work by Leclair and Mussardo [5] in the bulk case. 

In order to provide evidence for the conjecture, the proposed formula is developed in 
a low-temperature series, with the details described in appendix A. The low-temperature 
expansion of (|1.2p is then evaluated using an independent method developed in [6]. This 
approach requires the knowledge of boundary form factors in finite volume (up to corrections 
that decay exponentially with the volume). Section 4 presents the relevant results from the 
paper |^, and appendix B provides some further details on the evaluation of diagonal matrix 
elements. The calculation itself is presented in section 5, with a particularly complicated part 
relegated to appendix C. Section 6 is devoted to the conclusions. 

2 The boundary form factor bootstrap 

The relations satisfied by the form factors of a local boundary operator were derived in [1]. 
Compared to the equations in [3], the ones presented here are slightly generalized to allow for 
more than one particle species. Such an extension was first given in [8]; the derivation of these 
equations is straightforward using the methods of [4]. 

Here the equations are listed without much further explanation. Take an integrable bound- 
ary quantum field theory in the (infinite volume) domain x < 0, with N scalar particles of 
masses rrta {a = \ . . .H). As usual in two-dimensional field theory, asymptotic particles are 
labeled with their rapidities ^, and their energy and momentum reads 



Both the bulk and boundary scattering are assumed to be diagonal and given by the two- 
particle S matrices 

Sa^a.iOl - 02) = e'^'^^-^^'^-'^^ (2.1) 

(where Saia2{Gi — 62) are the two-particle phase-shifts) and the one-particle reflection factors 

Ra{e) 

satisfying the boundary reflection factor bootstrap conditions of Ghoshal and Zamolodchikov 
[1]. For a local operator 0{t) localized at the boundary (located at x = 0, and parametrized 
by the time coordinate t) the form factors are defined as 

a;...a^(^l)--- )^m|C(*)l^li--- i^n)ai...a„ = 

rpO //)' a' ■ a a N„-«mi(X;coshei-J]coshe ) 

-^<...a;;ai...a„(,^l'---'^m>^l>---)^n)e '^ ^ ^' 



using the asymptotic states introduced in [9]. They can be extended analytically to complex 
values of the rapidity variables. With the help of the crossing relations derived in [4j all form 
factors can be expressed in terms of the elementary form factors 



(010(0)1^1,..., e„),„ = F,%„ (01,..., 0„) (2.2) 



which can be shown to satisfy the following equations: 
I. Permutation: 



Ki...aiai+i...a„((^l^--- , ^i, (^i+1, ■ ■ ■ ,(^n) - (2-3) 

Saia,+ iiOi - ^i+l)-^ai...ai+iai...a„(^l'---.^i+li^j'--- ) ^n) 

II. Reflection: 

^ai...a„(^l.--- ,(^n-l,dn) = Ra„{Gn)Fa^ ^^{Oi, . . . ,9n-l,-9n) (2-4) 

III. Crossing reflection: 

F^^,,,,^^{01,92,...,en)=RaAi^-Ol)FZ..aJ'^i^-Ol,O2,...,On) (2.5) 

IV. Kinematical singularity 

-i Res FZa,...a„ {e + t7r,9',9,,..., 9^) = (2.6) 

Caa'U-flSaaA0-0^)SaaA9 + 9AFZ..aM,...,9n) 

where Caa' = ^aa' is the charge conjugation matrix (a denotes the antiparticle of species a). 

V. Boundary kinematical singularity 

n 



i^^^FZ,...aS0 + -,9r,...,9n) = '-^{l-\[SaaA^-9i))FZ...aSe,,...,9n) (2.7) 



where Qa is the one-particle couphng to the boundary 

There are also further equations corresponding to the bulk and boundary bootstrap struc- 
ture (i.e. bound state singularities of the scattering amplitudes S and R), but they are not 
needed in the sequel. The equations are supplemented by the assumption of maximum ana- 
lyticity i.e. that the form factors only have the minimal singularity structure consistent with 
the bootstrap equations. We remark that it is a general property of non-trivially interacting 
diagonal factorized scattering theories that their amplitudes are fermionic: 

SaaiO) = -1 

and as a result of eqn. (|2.3I) all form factor functions satisfy an exclusion property (Pauli 
principle), i.e. they vanish when any two of their rapidity arguments coincide, together with 
the corresponding species indices. 

It was shown in [10] that the space of solutions of the above equations is consistent with 
the operator spectrum predicted by boundary conformal field theory in the Lee- Yang and 
sinh-Gordon model. More recently the author gave a general procedure to construct solutions 
with a specific scaling dimension starting from an appropriate solution of the bulk form factor 
axioms |11| . 

We remark that using the bulk form factor bootstrap (cf. [12] for a review) as a guide it 
is straightforward to extend these axioms for non-diagonal scattering, i.e. particles with an 
internal degree of freedom. Some results for such theories (albeit only for diagonal boundary 
scattering) can be found in [131 [H] . 

3 A conjecture for the expectation values 

Consider a theory with a spectrum that contains a single massive particle species of mass 
m. Leclair and Mussardo proposed the following expression for the bulk finite temperature 
one-point functions |S|: 

(^>'' =t^,fi(r §T^^l /-w «.) (3.1) 



n! 

n=0 4=1 



where /!„ is the connected diagonal form factor of the local bulk operator A, R = 1/T in 
terms of the temperature T, and e{6) is the pseudo-energy function, which is the solution of 
the thermodynamic Bethe Ansatz (TBA) equation 



/OO ^Q> 
—(p{e-e')\og{l + e-<'^''^) (3.2) 



d9' 
e[tl) = mK cosh t) — j 

J —OO 

where 



is the derivative of the two-particle phase-shift introduced in (12.1(1 . The factor 1/n! takes 
into account the fact that a complete set of n-particle in-states is obtained with the ordering 



^1 > ^2 > • • • > ^n, but the integrals can be extended to the entire space using the fact that 
the functions /2n(^i> ■•■) ^n) are symmetric in all of their arguments. 

The main idea behind the formula (I3.ip comes from the TEA expression of the free energy 



/OO JQ 
— mcosh(0)log(l + e^ 



<eh 



which shows that the finite temperature vacuum can be considered as a free Fermi gas of 
quasi-particles for which the thermal weight is given by the pseudo-energy function e{0). The 
essential condition necessary for the validity of this picture is that the complete set of states 
used to derive (13. ip must be inserted at a position which is asymptotically far from any local 
operator insertion, so that their distribution is governed by the unperturbed finite temperature 
ground state. This is the reason why the Leclair-Mussardo conjecture does not work for the 
two-point functions ^5], because the states inserted between the two local operators cannot 
be asymptotically far from the positions of the operators which are themselves located at a 
finite distance from each other. 

From figure 11.11 it is obvious that a complete set of asymptotic states can be inserted at 
X = — OO where their distribution is unaffected by the presence of the boundary operator O. 
The only difference to the bulk case is that the complete system of in-states is spanned by 
multi-particle states with all their rapidities positive (i.e. with all particles moving towards 
the boundary), so the natural generalization of (13.11) is 



n=0 1=1 \ / 



where F^^ is the connected part of the diagonal form factor of the local boundary operator O: 

FIJ^l, ..., On) = (01, 02, ... , en\0{t = O)|0i, 02, . . . , 0„)-™-*-'^ (3.4) 



which is again symmetric in all their variables as a result of equation (I2.3|l . The precise 
definition of the connected matrix element (valid both for bulk and the boundary operators) 
is specified later in subsection 4.2. 

The conjectured expression (|3.3p can be checked against a calculation of the low-temperature 
expansion using the boundary form factors; this calculation is performed in the sequel. How- 
ever, the kinematical residue equation (|2.6p implies that diagonal matrix elements contain 
disconnected terms which are infinite, and therefore must be regularized. As shown in [6] a 
natural regularization can be obtained by putting the system in a finite volume, which was 
implemented for the bulk case in [l6l[6j and for the boundary case in ^. 

For completeness we note that the conjecture (13. 3p can be extended to a theory with 
multiple particle species and diagonal scattering in the following form: 



where 



i^a,...aj^l,-,^n)= ai...a„ (01,02,..., 0„|O(t = 0)101,02, ^ ^<^onnected 



ai...a„V"i, ■■■■,"n) — ai...an\^ V,"2, ■ ■ ■ ,^n\^\'^ — ^J\"L,^2, ■ ■ ■ ,i^n/ai...a„ 

while the pseudo-energy functions satisfy 

— V?a6(0-0')log(l + e-^«^(^')) 



where 

VabiO) = ^Sat{9) (3.6) 

are the derivatives of the two-particle phase-shifts introduced in (|2.ip . For the sake of simphc- 
ity the species labels will be omitted from now on, i.e. every formula will be written for the 
case of a single particle species; the extension to multiple species (with diagonal scattering) is 
rather straightforward. 

Eqn. (|3.3|) can be expanded systematically order by order in e""^^ which yields a low 
temperature expansion, following the procedure implemented for the Leclair-Mussardo formula 
(|3.1|) in [6]. The detailed calculation is performed in Appendix A with the following result: 

{Of = ai + a2 + a3 + (e"^™^) (3.7) 

where 

/ d8i /mi?, cosh Si „-2mi? cosh 6*1 , „-3mi?,cosh6'i\ rpc/n \ 

+\ r^ r^e-'"^(-^^^^+-^'^^^)ci>,2(F|(^i)+F|(^2)) 

1 f dOi f d92 f rf^3 -mi?,(cosh6li+cosh6l2+cosh6l3)/'^ ^ 

^2], ^k ^/o ^' ^^"^^^ 

+$12$23 + ^13^23) i^K^'s) 

/■°° dOi r°° d^2 mR(2coshei+cosh92) fojpc(a \ , ^ rpcfn \\ fl, 

2 7o 27r Jo 2tt 

-Jo -^L ^' "^^^'^''^^ 

+\ r^ rp r^e-^(™^^^^^+-^^^^+-^^^3)(ci>^2+ci,,3)i.c(^2,^3) 

2 Jo ^■^ Jo 2^ Jo ^■^ 

0-3 = I r°° dOi r°° de2 f°° de3^-mRicosh9^+coshe2+coshe;)pC(^0^^Q^^Q^^ 

6 Jo 27r io 27r Jo 27r 



and 



$ij = V7(^i-%) + <^(^i + ^j) 



For later convenience some terms were reordered by reshuffling the integral variables. 

In the sequel this result is compared to the result obtained from explicit evaluation of the 
finite temperature Gibbs average. In order to perform this calculation it is necessary to use 
finite volume as a regulator, and so now we turn to the issue of boundary form factors in finite 
volume, based on the results of [7]. 

4 Boundary form factors in finite volume 

4.1 Bethe-Yang equations 

Let us consider an integrable boundary quantum field theory with particles of species a = 
1, . . . , N and corresponding masses rua in finite volume L as shown in figure H?T1 As in section 



a 13 

I 1 I 

u y \j 

X = —L X = 

Figure 4.1: The setting of Fig. Il.ll in finite volume 

2, the bulk and boundary scattering is assumed to be diagonal and given by the two-particle 
S matrices 

and the one-particle reflection factors 

i?(")(0) = e*^«"^W , 4/^)(0) = e^'5^"W (4.1) 

where a and 13 denote the left and right boundary conditions, respectively. 

In the diagonal case, the multi-particle energy levels in a finite volume L are described by 
the following Bethe-Yang equations [TT] : 

Q.(^l>--->^n)ai...a„ = 2^^' (4-2) 

where the phases describing the wave function monodromies are 

Q.(^i,---,^nU...,„ = 2m„^,Lsinhe,+4^)(0,)+5(f (0,) 

Here all rapidities 6j (and accordingly all quantum numbers Ij) are taken to be positivqj. The 
corresponding multi-particle state is denoted by 

li-'l! • • • ) ^n\)ai...an,L 

and its energy (relative to the ground state) is 



where \ 6j \ is the solution of eqns. (|4.2|) in volume L. The energy calculated from the 

L J j=l,...,n 

Bethe-Yang equations is exact to all order in 1/L; only finite size effects decaying exponentially 
with L are neglected. 

^Boundary reflections change the sign of the momentum, so flnite volume multi-particle states can be 
characterized by the absolute value of the rapidities. 



4.2 Matrix elements in finite volume 

In general infinite volume and finite volume matrix elements are just related by the square root 
of the ratio of normalization of the corresponding states [Jl [16]. This results in the following 
relation: 

b,...b^{{l[, ■■■, 4}|O(0)|{/i, . . . , In})a,...a^,L = 

br^- b^-^--^\"' ^ ' ^ + ©(e-A*^) (4.3) 

yPai...a„ (Ol,---, On)Pbi...bm 0'i, ■ ■ ■ , O'm) 

where F^^ ^^{6i, . . . ,6n) is the form factor of the operator O (in the infinite volume theory, 
i.e. on the half-line x < Q)a9j\ is the solution of eqns. (14. 2p in volume L for the set 

L ) j=l,...,n 

of quantum numbers {Ii, . . . , /„} (similarly for < 6'- > and {/{, . . . , ImDi ^^d 

la a \ A i- I ^^k\yi^ ■ ■ ■ , "n)ai...an I /4 4\ 

Pai...a„(6'i,...,6l„) = det<^ — ^ (4.4) 

is the finite volume density of states, which is the Jacobi determinant of the mapping between 
the space of quantum numbers and the space of rapidities specified by the Bethe-Yang equa- 
tions (14.21) . An explicit expression for the derivative matrix of the Bethe-Yang equations (|4.2p 

is 

^ = -^a,a,{e,-eu) + ^a,a,{ej+eu) , j^k (4.5) 

where 

are the derivatives of the boundary phase-shifts defined in (|4.ip . while the if are the derivatives 
of the bulk ones as written in (13. 6|) . 

Eqn. (|4.3p is valid as long as the sets of the rapidities corresponding to the two states, 

\ Oj \ and \ O'j \ , are disjoint i.e. when there are no disconnected contributions. 

L J j=l,...,n L J j=l,...,m 

For diagonal matrix elements 

a^...a„{{h, ■ ■ .,In}\Om{h,. . . , In})a,...a„,L 

a more careful analysis is required [6l [7]. According to (14. 3p for this case it is necessary to 
consider 

-Fa„...aiai...a„(6'n +i7r,...,9l + ITT, 9i, ..., 0n) 

Because of the kinematical poles the above expression is not well-defined. The bulk kinematical 
singularity axiom (12.61) implies that the regularized version 

Fa„...aiai...a„{On + in + e„, ...,0i + ilT + ei, ^i, ..., 0n) 



has a finite limit when ej — > simultaneously. However, the end result depends on the direction 
of the limit, i.e. on the ratio of the ej parameters. The terms that are relevant in this limit 
can be written in the following general form: 

Fa„...aiai...a„i6n + i-n- + €n, ■■.,01 + ITT + ei, 9i, ..., 0n) = (4.6) 

n ^ n n 

• 1 ^* • 1 -1 

where A^^"f^ is a tensor of rank n in the indices ii,... ,in which is symmetric under the 
exchange of indices that correspond to particles of the same species, and the ellipsis denote 
terms that vanish when taking ej — > simultaneously. 

The connected matrix element can be defined as the ej independent part of eqn. (14. 6|) . i.e. 
the part which does not diverge whenever any of the ej is taken to zero: 

F^^,„,jeu:.,Bn)= Y. -^m:::p:(^i'---'^") (4-7) 

(Pl---Pn) 

where the summation goes over all permutations (pi, . . . ,pn) of the numbers l,...,n. As 
shown in appendix B, all other evaluations of the diagonal matrix elements (|4.6|) can be 
readily expressed in terms of the connected amplitudes. 

It was shown in [7j that a natural generalization of an expression proposed earlier by Saleur 
[T5] for bulk diagonal matrix elements can be extended to the boundary case in the following 
waja: 

a,...a„({/l . . . In}\Om{h . . . /„})a,...a„,L = (4-8) 

^ TT E KiA)imkeA)Pa,...aAOu ■ ■ ■ A\A) + 0(e-^^) 

Pai...a„it^l,---,t^nj Ac{l,2,...n} 

The summation runs over all subsets A of {1, 2, . . . n} and again \ 9j > is the solution 

I ) j=l,...,n 

of eqns. (|4.2p in volume L for the set of quantum numbers {/i, . . . , /„}. For any such subset 
the corresponding species index list is defined as 

a{A) = {ak}keA 
and 

Pai...a„(^l, . . .,On\A) = det J^'-"" {01, . . . , 0„) (4.9) 

is the appropriate sub-determinant of the n x n Bethe-Yang Jacobi matrix 

T ff) o \ _ (^QkyJl, ■ ■ ■ ,"n)ai...a„ ,. ^ „x 

Jai...aA9li---,Bn)kl — ^^ (,4.iUj 

obtained by deleting the rows and columns corresponding to the subset of indices A. The 
determinant of the empty sub-matrix (i.e. when A = {l,2,...n}) is defined to equal 1 by 
convention. It is also shown in appendix B that the symmetric evaluation which gave a very 
convenient alternative to ()4.8p in the bulk ^, behaves rather differently in the boundary case. 



^Note that here the original result of [7] is extended to the case of several particle species. 



5 Expansion of finite temperature expectation values 

5.1 Low-temperature expansion for one-point functions 

The procedure leading to a well-defined low-temperature expansion was outlined in section 
7 of [6]; details about the validity of the method and the existence of the limits taken are 
omitted (the interested reader is referred to the above paper for details). Let us evaluate the 
finite temperature expectation value of an operator O located at x = in a finite but large 
volume L, according to the setting introduced in section 4: 



TrL {e-^^L 



(Q)l= r^\.-RH.. ,T = l/R (5.1) 



Hl is the finite volume Hamiltonian, and Ttl means that the trace is now taken over the finite 
volume Hilbert space. The expectation value (O) can be recovered in the limit L — > oo which 
means that the left boundary condition a in figure 14.11 plays an auxiliary role, and the end 
result can only depend on the x = boundary condition /?; this issue will be taken up again 
in subsection 5.5. 

In the calculation below particle species labels are dropped for simplicity (they can be 
easily reinstated if necessary) and we use the simplified notation F2n for the n-particle diagonal 
matrix element introduced in (|3.4p . It is also convenient to introduce a new notation: 

1^1, • • • , On)L = K-^l, • • • , In})L 



where 9i, . . . ,9n solve the Bethe-Yang equations (|4.2p for n particles with quantum numbers 
Ii, . . . , In in volume L; as remarked in subsection 2.1, all of the rapidities can be taken positive. 
The low temperature expansion of (|5.1|) can be developed in orders of e""*^ using 

TVL(e-^^^O) = (O)i + ^e-'"^™'^^'''''(0«|O|0W)L 

g(2) g(2) 



+ 1 ^ '^-™ii(cosh<)+cosh4^)+co.hef))^^{3)^^(3)^^(3)|^|^(3)^^(3)^^(3)^^ 

/)(3) fl(3) „(3) 

+0(e-^"^) (5.2) 



and 



2 

9(1) g(2)^g(2) 



f, Z — ^ 



6 

/,(3) n(3) fl(3) 
"l .''2 't^3 



10 



The denominator of (15. 1|) can then be easily expanded: 



TVL(e-«^i; 




e 



V 

I _ 1 V^ ' -m/?(cosh6»f^+coshef^) 




-mR cosh 61(1) | , / ST^ -mRcoshOW \ V^ ' -mH(cosh6lf^+cosh0^^') 



_l_ I V^-mR cosh 61(1) 1 V^ '^ 

V(i) / ef),ef 

1 V^ ' -m_R(cosh6lf +cosh6»^^'+cosh6»*^') _|_ ^/ -4m_R\ /r a\ 



q(3) fl(3) ^(3) 



The primes in the multi-particle sums serve as a reminder that there exist only states for which 
all quantum numbers are distinct. Since it was assumed that there is a single particle species, 
this means that terms in which any two of the rapidities coincide are excluded. All n-particle 
terms in (|5.2p and (|5.3|) have a 1/n! prefactor which takes into account that different ordering 
of the same rapidities give the same state; as the expansion contains only diagonal matrix 
elements, phases resulting from reordering the particles cancel. It is also crucial to remember 
that in the boundary case the summations only run over positive values of the rapidities (cf. 
section 3). The upper indices of the rapidity variables indicate the number of particles in the 
original finite volume states which helps to keep track which multi-particle state density is 
relevant. 

It is also necessary to extend the finite volume matrix elements to rapidities that are 
not necessarily solutions of the appropriate Bethe-Yang equations. The required analytic 
continuation can be written down using eqn. 



^"(^^'•••'^")^ Ac{t^...n} 

(5.5) 

where the volume dependence of the n-particle density factors was made explicit and the form 
factors are computed from solutions of the bootstrap equations in section 2 with the boundary 
condition /3. It is apparent that the continuation is specified only up to terms decaying 
exponentially with the volume L but this is sufficient for the evaluation of the L — > oo limit 
ofKTh. 

It is useful to notice that unitarity and real analyticity imply that all the phase-shift 
derivatives 

^{0) = ^{9) , tl;^"\9) = ^5(")(0) , i.(^\9) = ^6^P\9) (5.6) 

are real and even functions. Another important observation is that the exclusion principle 
(cf. section 2) implies that the amplitudes F|„(6'i, . . . ,9n) vanish whenever any two of their 
rapidity arguments coincide. In addition, the connected form factor functions are symmetric 
under permutations of their arguments according to their definition (|4.7p and (in contrast to 
the bulk case) they are even functions in all of their rapidity arguments separately i.e. 

-^2n(^l;^2, ■ ■ ■ ,9n) = -^2n(~^l; ^2, • ■ ■ ,9n) 

which is a result of the reflection equation (|2.4|) satisfied by the form factors. 

11 



5.2 Lowest order terms 

The leading correction is 

{o)1 = {o)l + Y, 6-"^^"°^^^"' [{e^^\o\e^^^)L - (o)l) + oie-^"^^) 
0(1) 

Prom JSSl) 

{e\0\9)L - {O) = ^fm + O (e-^^) 

Note also that the difference between the finite volume vacuum expectation value and the 
infinite volume one decays exponentially with L 

{0)l -{0)^0 (e-^^) 

Prom now on such exponential corrections will simply be omitted. In the large L limit the 
summation can be replaced by the integral 



and therefore 



E^/|^..(«) 



'0 

5.3 Corrections of order e"^"^^ 
To this order one has 



OO JO 

(0)« = {0)+ I i!^i7|(0)e-'"«™«he ^ o(e-2™«) (5.7) 

27r 



em 



'1 / \"2 



I 



g(2) g(2) 

Using the symmetry of the first term in the rapidities and separating the diagonal contribution 
from the double summation on the last line leads to 



2 

g(l) 



9(1) 



+ 



1 ^ e— ^(™^t^<^+™^^^^2'')((0f\^f|O|0f\0f)i-(O)i 

g{2) g{2) 

1 ^ e-2-^™^^^^^^((^f\^f)|O|0f\^f)),.-(O)z.) 



9(2) ^g(2) 



12 



The terms containing two independent rapidity sums can be written as 

+F4(«i,''2) + ««i,«2|{l})F|(«i)+p(9i,''2|{2})F|(92)j 

where 

pi [6) = 2mL cosh + V^"^ (0) + V'^'^^ (^) 

is the one-particle state density, while 

p(^i,^2|{l}) = 2mLcoshe2 + iJ^''\92) + i^^'^\02)+vi92-0i)+^{92 + ei) (5.8) 

are the corresponding sub-determinants of the two-particle Bethe-Yang Jacobian, evaluated 
according to (|4.9p . Taking L — > oo 



5^f = 11 |i^e-™^(^°^^^i+=°^^^^) [F|(^i, ^2) + (^(^1 - 92) + ^{9, + 92)) {F^{9,) + ^2^(02))] 
The diagonal contribution contains a single rapidity sum 

g{2)^g{2) 

for which one needs to evaluate the density of states for a degenerate two-particle state. The 
appropriate Bethe-Yang equation reads 



2mLsinh0f ^ + 5{0) + S (29'i^j + J^") (^^f ^j + J^'^) (^^f ^j = 27r/i (5.9) 

and so the summation can be replaced by 

pi2{9) = 2mLcosh9 + 2ip{29)+i:'-''\9i)+4^^^\9i) (5.10) 

On the other hand from (15.511 it follows that 



{9,9\0\9,9)l-{0)l = -^^^Fi{9u92)+pi0i,02\{l})F^{ei) 

P2{9,t)) I 

+p{9u02\{2})F^{92)]^ ^ +0{e->^^) 

where the p are given in (I5.8|l and the factor at the front can be calculated from (|4.4|) 

P2{9, 9) = Am^L^ cosh^ 9 + 0{L) 
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In addition, the exclusion property can be used to substitute F^{9,9) = 0. Taking the hmit 
L ^> oo results in 

and so the total contribution at this order reads 

oo ja 

"'^ „-2mi?,coshe J7C//)N 

^e F,{6) 

1 /■ dOi^ r d92 rrai?(coshei+cosh6>2) 



2 7q 27r Jq 2tt 



Fi{ei,e2) 

+ iviOi - 02) + ^{01 + 62)) (F|(0i) + F|(02)) 1 (5.11) 



5.4 Corrections of order e 



-3mR 



L 



dOi f d6'2^m^(2 cosh ei+cosh 6(2) 



27r Jo 27r 



This calculation proceeds in a similar way but it is rather long and so it is relegated to appendix 
C. The net result is 

6 Jo 27r Jq 2-k Jo 2tt 

+3F|(6'2, ^3)(^12 + ^13) + 3F|(03)($12$13 + $12^23 + ^13^23)] 

"i^|(^l,^2) 

2F|(0i) + -F^{e2)\ $12] + j ^e-3-«^°^^^^F|(^i) (5.12) 

where ^ij = ^piOi - %) + v{9i + 0j). 

5.5 Discussion of the results 

It is very important to note that in the order by order corrections (|5.7|) . (jS.lljl and (|5.12|) . 
the dependence on the boundary condition /3 at j; = is only carried by the form factors 
Fj^. However, in the intermediate calculations the Bethe-Yang determinants enter, which 
depend on the boundary conditions a and /3 in a symmetrical way: according to eqn. (14.51) . 
the boundary phase-shift derivatives always appear in the combination 

^(")(0)+^(/3)(e) 

The fact that all such terms drop in the L — > oo limit is necessary for consistency since the 
end result can only depend on the boundary condition /3 imposed at x = 0, but not on the 
auxiliary (and indeed arbitrary) boundary condition a imposed at x = —L (cf. figure HH]). 
Summarizing the results, the expansion of the one-point function reads 

(O)^ = (O) + Si + S2 + S3 + O (e-^'"^) 
where 



r"00 

-mR cosh 9 
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while S2 and S3 are given in eqns. (15. lip and (|5.12l) . respectively. Note that this result 
exactly coincides with the expansion (13. 7p of the conjectured formula (|3.3p . which is a strong 
reason to believe that the conjecture is indeed correct to all orders (especially in view of the 
very nontrivial structure of the third-order correction terms). 

There is a rather obvious structural similarity between the bulk formula (jS.ip and the 
boundary one (|3.3p . Taking into account the symmetry of the pseudo-energy function exploited 
in appendix A, it is possible to bring the bulk and boundary cases into correspondence by 
interchanging the following ingredients: 

bulk boundary 



TOO dOi 

J— 00 2tt 


(■00 ddi 

h 2-K 


nnie,,...en) 


FUe,,...e^) 


^{e, - Ok) 


ipiOj - 9k) + v{9j + Ok) 



(some care must be taken on the second line to follow properly the particle labels of /|, since 
the bulk connected two-particle form factor is actually independent of the rapidity and thus the 
argument is usually omitted). Since Theorem 1 of appendix B is related to the corresponding 
bulk theorem of [6] via the correspondence implied by the last two lines in the above table, it 
is also possible to express the expansion (13. 7|) in terms of symmetric form factors analogously 
to the result obtained in [6]: 

—mi?, cosh —2mRcoshd , —ZmRcoshd 



/■oo jn 

{Of = {o)+j^ -Fm 



4r^r^^^'<-«^ 



-~vaR(coshdi+cosh92) q — mi?(2cosh9i+cosh62) 

1 /■ d^i /■ d^ /■ '^^ps/Q Q \g-miJ(cosh6li-|-cosh6l2+cosh6l3) 

6 io 27r 7o 27r 7o 27r '^ 

/■ dOi r 6,62 j^sfn \ ^T^s/n\ (f) „-mR(2coshei+coshe2) I n ^o-4mi?\ 



where eqns. (|B.5|B.6|B.7|) were used, together with the freedom to relabel some integration 
variables. However, note that this is not automatically guaranteed in the finite volume formal- 
ism used in the present section, since the computation makes use of the various Bethe-Yang 
determinants which depend explicitly on the combination ip{9j — 9^) — (p{9j + 9^) as pointed 
out in appendix B. The agreement between (|3.7|) and the corrections in eqns. (|5.7|5.11|5.12p 
shows that this dependence drops out after the limit L — > 00, which is far from trivial, albeit 
required for overall consistency. 

6 Conclusions and outlook 

The main result of this paper is eqn. (13. 3p (or its generalization (13. 5p ) which provides a way to 
evaluate finite temperature expectation values of boundary operators in terms of form factors. 
At first sight all the rest of the argument (i.e. the low-temperature expansion using the 
finite volume regularization) is only developed in order to verify this conjecture. However, as 
already pointed out for the bulk case discussed in [6], the finite volume regulator can be used 
to evaluate two-point (or even higher) correlation functions at finite temperature. There has 
been some development in the bulk case [3 [l5l [IBl [IH [20] , but there is a general problem that 
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the regulator imposed to deal with the disconnected contributions is rather ad hoc. The failure 
of Delfino's proposal for the bulk finite temperature expectation values [221 EB] shows that 
the ambiguity inherent in the regularization procedure (which is manifested in the directional 
dependence of the diagonal limit discussed in subsection 4.2 and appendix B) must be taken 
seriously. 

However, as pointed out already in [6j, finite volume as a regulator is guaranteed to give 
a correct answer as a matter of principle, since it provides a physical way to regularize the 
form factors entering the expansion. Therefore it would be very interesting to apply the ideas 
presented in J6] and here to compute bulk and boundary two-point functions, respectively. 

Another interesting issue is to obtain an extension of the finite volume description of form 
factors to non-diagonal scattering theories, both in the bulk and on the boundary. Since the 
description of finite volume energy levels is known and is not very complicated (one obtains 
scalar Bethe-Yang equations after suitably diagonalizing a family of commuting transfer ma- 
trices, cf. [21] and references therein), it can be expected that the necessary description of 
form factors is not too difficult to find. One can then use these results to evaluate finite 
temperature averages and correlators in the non-diagonal case as well. 
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A Low-temperature expansion of the conjectured formula ([3TT]) 

First the pseudo-energy function e{6) must be expanded to the necessary order. Using the 
fact that e{9) is an even function, the TBA equation can be written in the form 



/•OO J of 

e{9) = mRcosh{9) - / t^ [f{e - 6') + ^{9 + 6')] log(l + e" 



<d'Y 



Iterating this equation twice with the starting value e^^>{9) = mRcosh{6) and taking care to 
expand the logarithm one obtains 

e(^i) = RE,- r^'^ue-^^^-l r^<^ue-'^'^ 
Jo ^TT ^ Jo ^■^ 

Ei = mcosh9i , ^ij = (f{9i - Oj) + (f{9i + 9j) 



where 

which leads to 



Jo 27r 
' f^ d92^ 
Jo ^' 

+e-«^^ r p r ^^^^^^^^-RiE.+E,) ^ o ^^-AmR^ (A^^ 

Jo 2^ Jo 2^ 



10 

OO J/)_ /-OO 
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Recall that (13. 3|) reads 



n=0 i=l \ / 

Using (lA.lh and the geometric series 

— = e"' - e"^' + e"^^ + . . . 

1 + e-^ 

this can be expanded in orders of e"*"^. One obtains 

(C))« = ai + a2 + a3 + (e"^'"^) (A.2) 

where 

a, = r^ (e-«^^ - e"^«^^ + e'^^^^^) F^{9,) + T ^ T ^e-^(^^+^^^^^,F^{e,) 
Jo ^^ Jo ^^ Jo ^^ 

^ Jo ^^ Jo ^■^ Jo ^■^ 

2 70 27r io 27r ^^ '' '^ 

io 27r 7o ^vr 

are the one/two/three-particle contributions expanded to 0{e~^"^^). 

B Relation between different evaluations of the diagonal matrix 
element 

Here the arguments of [6] are generalized to the case of boundary form factors. The goal is to 
compute the general expression 

-^ai...a„(6'l,... ,9n\f-l,...,en) = i^a„...aiai...a„ (^'n + ivT + e„, ..., 6*1 + ITT + €i, 9i, ..., 9n) (B.l) 

for infinitesimal values of the e^. It is also interesting to consider the symmetric evaluation 

K...aA^l, ■■■,0n) = lim Fa„...aiai...aJ^n + iTT + e, ..., 9i + iTT + €, 9i, ..., 9n) (B.2) 

Let us take n vertices labeled by the numbers 1,2, ... ,n and let G be the set of the directed 
graphs Gi with the following properties: 

• Gi is tree-like. 

• For each vertex there is at most one outgoing edge. 
For an edge going from i to j we use the notation Eij. 

17 



Theorem 1 (IB.ll) can be evaluated as a sum over all graphs in G, where the contribution 
of a graph Gi is given by the following two rules: 

• Let Ai = {«!, 02, • • • 1 Om} be the set of vertices from which there are no outgoing edges 
in Gi. The form factor associated to Gi is 

Ko,^ ...a^^ (^ai , 6'a2 , • • • , ^a„ ) (B.3) 

• For each edge Ej^ the form factor above has to be multiplied by 



'' ^ik 



^cl>. 



efc 
where 

^jk = "Paja^iOj - Ok) + Va,ak{Oj + Ok) = ^kj 

Proof The proof goes by induction in n. For re = 1 there is only a single way to take the 
limit and so 

Fa{Ol\ei) = F^iOi) = F-aa{i7T + 0^0^) 

This is in accordance with the theorem, because for n = 1 there is only the trivial graph which 
contains no edges and a single node. 

Now assume that the theorem is true for n — 1 and let us take the case of re particles. 
Consider the residue of the matrix element (IB.ip at e^ = while keeping all the £« finite 

R = Res Fa^,„a„{Oi..On\ei..en) 

e„=0 

According to the theorem the graphs contributing to this residue are exactly those for which 
the vertex n has an outgoing edge and no incoming edges. Let Rj be sum of the diagrams 
where the outgoing edge is Enj for some j = l,...,re — 1, and so 



n-l 

J 






The form factors appearing in Rj do not depend on On- Therefore one gets exactly the diagrams 
that are needed to evaluate F2{n-i)(^i--^n-i|ei--en-i), apart from the proportionality factor 
associated to the link E^j and so 

-^i ~ ^j^jnFai...a„^iiOi..On-l\^l--^n-l) 

and summing over j yields 

R = (ei$i„ + e2$2n + • • • + en-l<^n-ln)Fai...a„-i{Ol--On-l\ei..en-l) (B.4) 

In order to prove the theorem, one only needs to show that the residue indeed takes this form. 
On the other hand, using the kinematical residue axiom (12. 6|) 

(n-l 
1 - JJ Sa„aj {On " Oj)Sa„aj {On " Oj -ITT- ej)Sa^a, {On + Oi)Sa„a, {On + Oj + ZTT + Cj 
i=l 
X-^ai...a„_i(6'l..6'„_i|ei..e„_i) 



• • •-^ • • ^-» 

12 12 12 

(a) (b) (c) 

Figure B.l: The graphs relevant for n = 2 



12 3 12 3 

(a) (b) 



12 3 12 3 

(c) (d) 

Figure B.2: The graphs relevant for n = 3 

which is exactly the same as eqn. (IB.4J) when expanded to first order in ej. 

Therefore the procedure described in the theorem gives the correct result for the terms 
that include a l/e„ singularity. Using symmetry in the rapidity variables this is true for all the 
terms that include at least one l/e^ for an arbitrary i. There is only one diagram that cannot 
be generated by the inductive procedure, namely the empty graph. However, there are no 
singularities (l/cj factors) associated to it, and it is identical to F2„(0i, . . . ,6n) by definition. 
Qed. 

Let us now illustrate how the theorem works for the case of a theory with a single particle 
species. In this case one can use the simplified notation introduced in (13. 4|) and similarly 
denote 

F2„(6'i,...,6'„|ei,...,e„) = F(6'„ + ivr + e„, ...,6'i + in + ei,9i, ...,0n) 

The case n = 1 is trivial: 

Fm=Fi(,e) (B.5) 

For n = 2, there are only three graphs, depicted in figure IbTTI Applying the rules yields 

Fi{9i,e2\eue2) = Fi{9i,e2) + $12 f-i^|(^2) + -F^{9 

which yields 

F|(^i,02) = Fi{e,,92) + $12 {F^{92) + Fi{9i)) (B.6) 

upon putting ei = 62. For n = 3 there are 4 different kinds of graphs, the representatives of 
which are shown in figure IB.2t all other graphs can be obtained by permuting the node labels 
1,2,3. The contributions of these graphs are 
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(a) : ^6^(01, ^2, ^3) 
(6) : ^-^<^i2Fi{92,e3) 

(c) : ^-l^-^^^^^^r^F^{e3) = ^-^^i2^23F^{e3) 

ei 62 ei 

(d) : !^f^cDi2$23i^|(^2) 

Adding up all the contributions and putting ei = €2 = £3: 

-^6(^1' ^2, ^3) = -^6(^1' ^2, ^3) 

+ (^12 + ^13)i^|(^2,^3) + ($12 + $23)^^4(^1,^3) + ($13 + $23)i^|(^l>^2) 
+ (F2^(ei)+F|(e2)+i^|(^3))($12$13 + $12$23 + $13$23) (B.7) 

It can be seen that these results are a natural generalization of the bulk ones obtained in |E] 
with if replaced by $. It is also important to keep in mind that contrary to the bulk situation 
F2 depends on the rapidity (in the bulk it is a constant since Lorentz invariance entails that 
all form factors depend only on rapidity differences). 

Now the finite volume diagonal matrix elements (|4.8p can also be re-expressed in terms of 
the symmetric evaluation. The first nontrivial case is n = 2 for which 



({/i,i2}|o(o)|{/i,/2}) = ,; -. iFi{eue2) + piOi,e2\{i})F^{ei) 



1 

P2[0l,f^2) 

+p{9^,e2\{2})F^{e2)) + (O) + 0(6-^^^) 



where ^1,^2 are the solutions of the 2-particle Bethe-Yang equations with quantum numbers 
h,l2- p denotes the appropriate sub-determinants (|4.10|) of the two-particle Jacobian matrix, 
while pn is the full n-particle Jacobi determinant (|4.4|) . It is straightforward to verify that 

Fi{ei,92) + m,O2\{i})F^{ei) + p{Ou02\mF^ie2) = 

F|(01,^2) +/Ol(^l)i^|(^l) +Pl(^2)F|(^2) 

+2ip{9i + 62) (FK^i) + Fi{92)) (B.8) 

The term on the last line shows that the analogue of Theorem 2 in [6] (which would make the 
expressions on the first two lines identical) fails in the boundary case. This results from the 
fact that the derivative matrix (|4.5p of the Bethe-Yang equations (|4.2|) carries a dependence 
not only on the combination ip{9j — 9^) + ^{9j + 9^), but also on ^{9j — 9^) — ^{9j + 9k). 

The relations (IB.6p . (|B.7|) and (IB.Sh were also verified numerically using the explicit form 
factor solutions presented in [7]. 
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C e ^"^^ corrections to the finite temperature one-point func- 
tion 

In order to keep the calculation manageable, let us introduce the following shortened notations: 

Ei = m cosh 9i 

{01, . . . ,en\0\9i, . . . ,en)L = {I ■ . .n\0\l . . .n)L 

p{ei,...,en) = p{i...n) 

p{9i,...,9n\{ai,...,ak}) = p{l. . .n\{ai, . . . ,ak}) 



Summations will be shortened to 



ei...9„ l...n 

E' - E' 



61... On l...n 

and for later convenience also denote 



^ij = ifiOi - Oj) + (fiOi + 9j) 



which satisfies $jj = $jj. 



Multiplying (|5.2p with (|5.4p and collecting the third order correction terms: 

}LY'e-RiEi+E2+E3) ((i23|C)|i23)i - {0)l) 
6 ■^— ' 

123 
V 1 J 23 



y. \ 1 / \ 2 / 12 J\3 / 

To keep trace of the state densities it is important to avoid combining rapidity sums. The 
constrained summations can be replaced by free sums with the diagonal contributions sub- 
tracted: 

E' = E-E 

12 12 1=2 

E' = E-(e + E + eV'^ 

123 123 \l=2,3 2=3,1 1=3,2/ 1=2=3 

where the diagonal contributions are labeled according to which diagonal the summation 
corresponds to, but otherwise the given sum is free, e.g. 

E 

1=2,3 
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shows a summation over all triplets 9[ ,62 ,9^ where 9[ = 62 and ^3 runs free (it can 
also be equal with the other two). Finally denote 

F{12...n) = F^^{9i,...,9n) 

so from (|5.5I) the necessary matrix elements can be written in the form 

p(123)((123|0|123)L-(0)i) = F(123) + p(123|{l, 2})F(12) 

+p(123|{l, 3})F(13) + p(123|{2, 3})F(23) 
+p(123|{l})F(l) +p(123|{2})F(2) +p(123|{3})F(3) 
p(122) ((122101122)^ - {0)l) = 2p(122|{l, 2})F(12) + p(122|{l})F(l) + 2p(122|{2})F(2) 
p(lll)((lll|0|lll)i-(0)i) = 3p(lll|{l})F(l) 

p{12) i{12\0\l2) l-{0)l) = F(12)+p(12|{l})F(l)+p(12|{2})F(2) 
p{U){{U\0\11)l-{0)l) = 2;5(11|{1})F(1) 

pil)i{l\0\l)L-{0)L) = F(l) (C.l) 

where the exclusion property was already used to eliminate form factors with equal rapidity 
arguments. 

One can now proceed by collecting terms according to the number of free rapidity variables. 
The terms containing threefold summation are 



^3'^ = ^E«"''^''''"'''^'''^((123|0|123)l-(0)l)-^^^((23|0|23U-(0)l) 

123 1 2,3 

Y 1 2 3 1,2 3 y 

Replacing the sums with integrals 

E - /It"!') 

1,2,3 -^ 

and using (IC.ll) 

S.f = I /'^^^e-^(^i+^2+^«)(F(123) + 3/5(123|{2,3})F(23)+3^(123|{3})F(3)) 
6 J 2tt 27r 2tt 

- \ I ^^^'''''''''^'''^^'^P^'^ (^('^) + 2p(23|{3})F(3)) 
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where some of the integration variables were reshuffled. The result is 



6 Jo ^ Jo ^io ^"^ 



^(3) ^ i_ I (m I 002 I a03 ^-m/?(cosh ei+cosh ga+cosh 6>3) 



'Q\'^l,"2,t>3) 



+3^1(^2,^3) (^12 + ^13) +3F|(03)($12^13 + ^12^23 + ^13$23)J (C.2) 

(3) 

(to derive the term on the second line note that the -F(3) terms in the integrand of S3 can 
be symmetrized in 61 and 62 without changing the value of the integral). 

It is also easy to deal with terms containing a single integral. The only term of this form 
is 

4^) = - Yl e-^(^^+^2+-^3) ((123|0|123)l - {0)l) 
1=2=3 

(3) (3) (3) 
When all rapidities 9^ ,62 ,0^ are equal, the three-particle Bethe-Yang equations reduce 

td 

2mLsinh^p) + 26 (20^^) + 5(") (^f )) + 5^^) (e['A = 2nh 



Therefore the relevant state density is 

P123 (0) = 2mL cosh 9 + 4(^(20) + ^^"^ (9) + ^j'-'^^ (9) 
and 

4" = 5/^e-««=.fe3(9l)({lll|0|lll>l-(0>l) 

= /|ie--.p(l)^™F(l)^-._^/^e— «.Ffm) ,C.3) 
where it was used that 

when L — > 00. 

The calculation of double integral terms is much more involved. The contributions con- 
taining two rapidity summations are 



y.(2) 
^3 



"M ^ + ^ + ^ I 6-^^^^+^^+^^) ((123|0|123)i - {0)l) 
yi=2,3 1=3,2 2=3, ly 

+^EE«"''^''''''''"''''^ ((23|0|23)^ - {0)l) 
1 2=3 

+^ E E e-«(^^+^^+^^) ((3|0|3)^ - {0)l) (C.4) 



2 

1=2 3 



The density of partially degenerate two-particle states was already computed in (IS.lOp . but the 
density of partially degenerate three-particle states is also needed. The relevant Bethe-Yang 



^Just as in l|5.9ll there are also contributions of the form 5(0), but these can be absorbed into a redefinition 
of II. 
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equations arcl 

2mLsmh9i + 6{ei-e2)+S{9i + e2) + S{2ei) + 6^'''^ {ei)+6^'^^{9i) = 2ttIi 
2mL sinh 62 + 26{92 - ^i) + 2.5(^2 + ^1) + S^"^ {O2) + 5^^^ (^2) = ^i^h 

where the first and the third particles are put as degenerate (i.e. I-^ = Ii). The density of 
these degenerate states is then 

pi3,2(12) = det ( "^'^ ''^^ ) (C.5) 

m = 2LEi + 99(^1 - 62) + ^{Oi + 62) + 2(^(2^1) + V^"^(^) + ^^^\0) 

r22 = 2LE2 + 2(^(^1 - 62) + 2^{9i + 62) + V'^"^ [0) + i^^^^ (6) 

ri2 = -^{Oi - 62) + ifiOi + 62) , rsi = -2v9(0i - 02) + 2(^(01 + 02) 

where it was used that if{9) = ip{—9). Using the above result and substituting integrals for 
the sums, eqn. (|C.4p can be rewritten in the form 



1 f d6id92 ji(2E,+Eo)Pl3,2{'i^'^) 



-e 



2p(112|{2,3})F(12) 



6 J 2tt 2tt p{ll2) 

+2p(112|{l})F(l) + p(112|{3})F(2) + . . . 

+ i/^f e-<-«..,.,lW3)^n3) 

where the ellipsis denote additional contributions that can be obtained by cyclical permutation 
of the indices 1, 2, 3 from those explicitly displayed inside the square bracket. These three sets 
of contributions can be shown to be equal to each other by relabeling the integration variables: 



1 f d9i C?6'2 Ji(2£;i+E2) Pl3,2(12) 



2p(112|{2,3})F(12) 



2 J 2tt 271 p{ll2) L 

+2p(112|{l})F(l)+/5(112|{3})F(2) 

The terms containing F{12) contribute 

J 27r 27r 



where it was used that 



/>l3,2(12)^,^^^ir^ o^^ _i , r^fT-l^ 



p(112) 



-p(112|{2,3})=l + 0(L- 



■^Just as in lj5.9|l there are also contributions of the form 5(0), but these can be absorbed into a redefinition 
of /l. 
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+ 



which results from (jC.SP and 

p(112|{2,3}) = 2mL cosh 6*1 + 2ip{0) + 2v9(20i) + +ip^°'\e) + iP^^\e) 
The terms containing -F(l) and F{2) combine to 

/ l^^e-.--. (-^„u,U,) ..(.)..(!)«) ni) 

A straightforward (albeit tedious) calculation leads to 

-^gf|^/^(112|{l})+p(2)p-i2(l) ^^^[|y^ = -2{^{ei - 62) + ^{9, + 02)) + O {L-') 

-^^f{^p(112|{3})+/^i2(l) = -^i9,-e2)-^i0i+e2)+o{L-') 

Note that the individual terms in these sums are proportional to L but their contributions 
drops out. For a more detailed discussion of such "anomalous" density contributions the reader 
is referred to [6]. 

The total contribution in the L — > 00 limit turns out to be just 

_i_^g-„^/^(coshel+2coshe.)(2i7|(^^) + -F^{92)){ip{e^ - 62) + ifiOi + 62)) (C.8) 
Summing up the contributions (IC.2J1 . (IC.3p . (|C.7p and (|C.8j) the end result is 

6 Jo 27r io 27r Jq 2tt 

+3F|(6'2,e3)($12 + ^13) + 3F|(03)($i2$l3 + $12^23 + $13^23)] 
pco JO roo jn i 

- / ^ / -^e-"^^(2coshei+coshe.)[^^c(^^^g^) ^ (2F|(0i) + -F|(e2))1>i2] 



oo 



"^1 _-3m_R cosh 6*1 pc 



7o 27r 
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